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Abstract 

We present photoionisation modelling of the intrinsic absorber in the bright quasar HS 1603+3820. We construct 
broad-band spectral energy distribution using optical/UV/X-ray observations from different instruments as an input 
to photoionisation calculations. Spectra from Keck telescope show extremely high ratio of Civ to Hi for the first 
absorber in system A, named Al. This value, together with high column density of Civ ion, puts strong constraints on 
photoionisation model. We use two photoionisation codes to derive hydrogen number density at the cloud illuminated 
surface. Estimating bolometric luminosity of HS 1603+3820, from typical formula for quasars, we calculate the 
distance to Al. Either for constant density cloud (modelled by cloudy), or stratified cloud (modelled by titan) we 
were able to find a single photoionisation solution, assuming solar abundances, which explains both the ionic column 
density of Civ, and high Civ to Hi ratio. The derived location is as close as 0. 1 pc, and situates an absorber even closer 
to the nucleus than the possible location of the Broad Line Region in this object. The upper limit for the distance is 
sensitive to the adopted covering factor and the carbon abundance. Photoionisation modelling always prefers dense 
clouds with the number density no = 10 10 , 10 12 cm" 3 , to explain observational constrains of HS 1603+3820 intrinsic 
absorption. This number density is of the same order as the number density in the disk atmosphere at the implied 
distance of the Al. Our results thus show that the disk wind escaping from outermost accretion disk atmosphere can 
build up dense absorber in quasars. 

Keywords: (galaxies:) quasars: absorption lines, atomic processes, radiative transfer, (galaxies:) quasars: individual 
(HS 1603+3820) 



1. Introduction 

High resolution optical/UV spectra show that over 20% of quasars exhibit broad absorption lines (BALs) from 
ionised species of heavy elements. In some objects lines are from high-ionisation species, as Civ, Nv, Ovi forming 
subclass of HiBALs, and in other cases absorption features originate from low-ionisation species as Mgn, Alii, Sin, 
Fen forming LoBALs subclass. 

In all BALs, lines can exhibit very complex blushifted profiles, indicating the presence of several absorbing sys- 
tems moving toward an observer with different velocities, for example in QSO 2359-1241 (Arav et al., 2001). Narrow 
line components indicate usually velocities of the order of a few hundred km s" 1 , while velocities of some broad lines 
can reach even 50,000 km s _1 (Weymann, 1995). 

Similar absorption in UV is observed in Seyfert galaxies (Crenshaw et al., 1999). Moreover, in 50% of those 
objects, the matter at higher ionisation state, so called "warm absorber", is detected in X-ray band, for example 
in NGC 3783 (Kaspi et al., 2001). Observed spectral features from Oviii, Mgxn, Sixiv, and even Fexxv indicate 
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temperature of an absorbing gas to be of the order of 10 6 - 10 7 K. Typical "warm" X-ray lines are narrower than those 
of BALs, and blushifted by velocities of the order of few hundreds kms~' up to 10 3 kms~' (for review, see Blustin 
et al., 2005). 

Many questions arise about the origin and physical conditions of those outflows. The first obvious is: does the 
UV and X-ray absorption occur in the same wind? There are two examples showing that the line detected in X-rays 
can be fitted exactly by the same velocity profile as the line detected in UV for two active galaxies NGC 5548 and 
NGC 3783 (Kaastra et al., 2002; Gabel et al., 2003). In general, to answer this question we have to collect multiband 
optical/UV/X-rays observations of quasars with resolution high enough to compare line profiles from different energy 
bands. However, presently working satellites as Chandra and XMM-Newton are not able to detect X-ray absorption in 
distant quasars, and very often X-ray spectrum can be modelled only as a power-law, like in case of Chandra detection 
of HS 1603+3820 (Dobrzycki et al., 2007, hereafter Paper I). 

Physical conditions in UV or X-ray absorbing clouds can be studied by photoionisation modelling. For that 
purpose, a broad-band spectral energy distribution (SED) of a considered source is needed. Nevertheless, many 
observations do not cover EUV band. For this reason in several recent papers on UV outflows, authors have used 
standard Mathews and Ferland (1987) composite quasar spectrum for photoionisation modelling. It is widely known, 
that the SED of incident photons determines ionisation and thermal structure of the absorbing plasma (Rozahska et al., 
2008; Chakravorty et al., 2009). In this paper, taking into account multi-wavelength observations of HS 1603+3820 
from MMT, Keck and Chandra telescopes we have constructed the SED of the quasar spectrum. This SED is used as 
an input for our advanced photoionisation modelling of the system using codes cloudy and titan. 

Photoionisation modelling provides us with a possibility to determine the distance of the cloud from the nucleus, 
only if the SED is dominated by soft disk component, as shown in Rozahska et al. (2008). Photoionised models 
depend predominantly on the ionisation parameter, i.e. a combination of the surface number density and a distance 
from the source. Models illuminated by only hard X-ray power-law are degenerated for the gas number density up 
to 10" cirT 3 . In thic case, we are not able to differentiate between matter located closer to the nucleus with higher 
density, and more diluted gas being farther away. The transmitted spectra through such two clouds look the same, since 
their ionisation and temperature structure are identical (Rozahska et al., 2008). In case of HS 1603+3820, where the 
SED is dominated by soft UV emission over X-rays, this degeneracy is broken due to the fact that the Comptonization 
becomes not important comparing to the free-free processes (Rozahska et al., 2008). Therefore, the hydrogen number 
density at the cloud surface can be estimated from the line ratios, and this density together with quasar luminosity, 
and ionisation parameter gives us a distance of an absorber from the source of continuum. 

Our method is complementary to other two methods used for distance derivation, based on different ways of 
determining the density number at an absorbing cloud surface. In the first one, it can be done from variability data, 
assuming that observed spectral changes are due to variation of the ionisation state (Krolik and Kriss, 1995). This 
method allowed us to estimate the number density of the warm absorber in Seyfert 1 galaxy NGC 4055 from ROSAT 
data to be of the order of n = 10 8 cirT 3 , implying the location of the gas at about R = 0.0024 pc (Nicastro et al., 1999). 
Nevertheless, in case of NGC 3783 two groups using this method got opposite results. Netzer et al. (2003) have found 
clouds at three different ionisation levels located on 3.2, 0.63, and 0.18 pc from the nucleus, while Krongold et al. 
(2005) found two clouds on 0.0029, and 0.0004-0.0008 pc, using the same 900 ks Chandra observations. Outflowing 
velocities of that absorbers were of the order of 750 km s . 

The second method uses high quality UV spectra of LoBALs, since number density can be estimated by comparing 
the ionic column densities of excited metastable states to their ground states of the most popular ions Hei, Cn, Sin, and 
Fen (Korista et al., 2008; Moe et al., 2009; Bautista et al., 2010). In case of three quasars presented in above papers, 
estimated number density is always of the order of 10 3 - 10 4 4 cirT 3 , indicating the gas location at about 3-6 kpc. The 
velocities of those absorbers are in the range of 1400-7500 km s -1 . 

In this paper we perform photoionisation modelling of an intrinsic absorption in quasar HS 1603+3820 with the 
aim to determine physical properties of an absorbing gas. This HiBAL type object is interesting since the ratio of Civ 
to Hi was found to be extremely high - above 20 (Paper I). The most interesting absorbing systems are outflowing 
with velocities reaching 3000 kms~'. By comparison of model parameters to those indicated from observations 
of HS 1603+3820, we are able to find a single photoionisation solution, which can explain observed ionic column 
densities of the first cloud in system A, reported in Paper I and marked as Al here. 

We use photoionisation modelling to determine uniquely both the hydrogen number density at the illuminated 
surface of the cloud, and the ionisation parameter. Those two quantities combined with the source luminosity allow us 
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to compute the distance from the central nucleus to the absorbing system. We present such calculations in this paper 
and compare results with previous estimations done by Misawa et al. (2005, 2007). 

This paper is organised as follows: in Sec. 2 we present the source, and in Sec. 3 we review its optical/UV/X-ray 
observations. Photoionisation models are presented in Section 4. The location of an absorber and connection to an 
accretion disk atmosphere is discussed in Sec. 5 and Sec. 6 respectively. We summarise our results in Sec. 7. 

2. The source 

Recently observed high redshift quasar HS 1603+3820 with z e m = 2.54 indicates extreme richness of absorption 
lines in its spectra (Dobrzycki et al., 1999; Misawa et al., 2003, 2005, Paper I). 

The source was first discovered in the Hamburg/CfA Bright Quasar Survey (Hagen et al., 1995; Dobrzycki et al., 
1996). With the brightness B = 15.9, it is among the top few brightest quasars known at such redshift. However, the 
most striking feature of HS 1603+3820 is its absorption spectrum, which indicates the presence of ~ 50 individual 
absorption systems, with 30+ of them having velocities higher than 10000 km s , when corrected for quasar redshift 
based on emission lines. Many - perhaps all - of these clouds are physically associated with the quasar, particularly 
because some lines show variability in the timescales of a few years (Misawa et al., 2007). 

High resolution observations from Subaru telescope with the High Dispersion Spectrograph, reported by Misawa 
et al. (2003, 2005) enabled to resolve many of those systems into a number of narrow components. The authors 
grouped the absorbers into several systems, designated with the letters from A to H. 

All of those systems, and two previously unknown, were indicated in the observations taken by MMT and Keck 
telescope and reported in Paper I. Furthermore, the authors have noticed that some clouds in system A have a high 
ratio of Civ to Hi column densities, reaching ~ 20. System A, and other systems spatially close to the quasar can 
be used as a probe for the intrinsic emission of the quasar itself, since the conditions in the systems are undoubtedly 
heavily influenced by the quasar flux (e.g. Crenshaw et al., 2003; Gabel et al., 2005; Scott et al., 2005, and references 
therein). 

Additionally, in Paper I, X-ray Chandra observations of HS 1603+3820 were reported, allowing to calculate the 
relative optical-to-X-ray slope, a ox — 1.49, which is typical for a quasar at this redshift (Bechtold et al., 2003). 

3. Observational constrains important for photoionisation modelling 

To explain the properties of an intrinsic absorption in HS 1603+3820 we need to construct illuminated continuum 
which affects absorbing gas. The observational data used in this paper are described in detail in Paper 1. In subsec- 
tions below, we shortly remind the reader values of fitted parameters for those observations, and summarise what is 
important for our modelling. 

3.1. Continuum shape and normalisation 

The optical/UV data spectrophotometry was taken from Scott et al. (2000), and it was obtained on July 4, 1995 
with Stewart Observatory Bok Telescope. Our continuum fit results with power-law photon index Tuv = 1-36 ±0.1. 
Observed optical flux at A res , = 1450 A (E rest = 8.55 x 10~ 3 keV) is EF E (B) = 3.20 x 10~ 12 erg s" 1 enr 2 as reported 
by Scott et al. (2000) (see their Tab. 3). In Fig. 1, blue points represent overall normalised flux in optical/UV domain. 

The X-ray data were taken with the Chandra X-ray Observatory on 2002 November 29. We performed the fit 
assuming the intrinsic quasar spectrum to be a power law, with fixed Galactic absorption towards HS 1603+3820 of 
A'h = 1.3 x 10 21 cm~ 2 . The best fit gave the photon index of Tx = 1.91 + 0.20 and normalisation at the rest frame 
for 1 keV of A = (2.43 ± 0.79) x 10~ 4 photons cm s 'keV . The X-ray part of observed and modelled spectrum is 
presented in the Fig. 1, red squares. 

HS 1603+3820 has been pointed at with IUE/LWP on 1995-08-19, but not detected. The observation gave the 
upper limit of 1.5 x 10~ 15 erg s _1 crrT 2 A -1 in the 2200-3100 A (observed) wavelength range. This limit is marked 
in green on Fig. 1. We note that this wavelength range covers the inferred position of Lyman-limit absorption from 
associated absorbers. It is therefore possible that this limit does not reflect low QSO's intrinsic luminosity in this 
range, although the individual absorbers do not appear to have sufficient column densities. 
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Figure 1: Normalised multi-band spectrum of HS 1603+3820 from MMT and Keck (blue points), IUE upper limit (green arrow), and Chandra 
(red squares) in the rest frame of quasar. Solid line represents the best shape of modelled spectrum used in calculations of photoionised intrinsic 
absorber presented in Sec 4. 



Table 1 : Ratios of ionic column densities of Civ and Hi observed in intrinsic absorption system A 1 , obtained for various assumptions about the 
covering factor of Hi cloud and fitting methods for system Al (z a bs = 2.41941 ± 0.00002). The name of the system and the No. Line ID are taken 
after our Paper I. 



Sys. No. 


log(iV Civ ) 


c f (N Civ ) 


b(N Civ ) 


log(Hi ) 


c/(Hi) 


b(Hi) 


N Civ IN m 


Al* 


14.97 ± 0.04 


0.36 


72 + 3 


13.65 + 0.03 


0.61 


62 ± 10 


20.37 ± 2.70 


Al 


14.97 + 0.04 


0.36 


72 + 3 


13.87 ±0.05 


1.00 


145 ± 12 


12.4 ±2.10 


Al 


14.97 ± 0.04 


0.36 


72 + 3 


14.11 ±0.05 


0.61 


132 ± 10 


7.1 ± 1.20 


Al 


14.97 + 0.04 


0.36 


72 + 3 


14.34 ±0.05 


0.36 


75 ±9 


4.2 ± 0.70 



8 global fit for all lines of system A from Paper I 



We adopt those values for all photoionisation models in Sec. 4, and for an estimation of bolometric luminosity 
in Sec. 5. We present the overall continuum in Fig. 1. Note that different parts of the spectrum are taken at different 
epochs therefore SED constructed by us can contain systematic errors. Nevertheless, for the time being, there is no 
other way to find continuum for this quasar. 

3.2. Optical/UV spectroscopy and parameters from spectral lines fitting 

HS 1603+3820 has been observed with two instruments. On 2001 August 19 we used the Echelle Spectrograph 
and Imager (Sheinis et al., 2002) on the Keck II telescope, with the combined exposure of 4800 sec. The 0.5 arcsec 
slit was used, yielding the resolution of 45-50 kms~'. 

Full procedure of spectral line fitting of optical/UV HS 1603+3820 spectrum is presented in Paper I. The line 
profile fits were performed using VPGUESS/VPFIT programmes 1 . We were able to fit simultaneously several lines 
from different ions, but characterised by the same redshift. 

In several clouds belonging to an intrinsic absorption systems the column density of Civ ion is bigger than the 
column density of ionised hydrogen Hi . The cloud Al was an extreme example, with the ratio of Na v IN Hi equal 
20 ± 3. In the present paper we concentrate on modelling this component in detail. 

In Paper I, the covering factor for the hydrogen and for carbon fitting were found to be different (0.36 and 0.61, 
correspondingly) which could argue against a single well defined cloud responsible for Al absorption. However, 
the kinematic width of the lines were similar. Therefore, we reconsidered the spectroscopic data by performing fits 
specifically for the component Al, instead of assuming the same covering factor for all clouds in component A. The 
old result and new results are given in Tab. 1 . 



'http://www.eso.org/~jliske/vpguess/, http://www.ast.cam.ac.uk/~rfc/vpfit.html 
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Figure 2: The fit of the hya region for Al component taking covering factor equal 0.36, as for Civ . Broad green dotted line represents an Hi line 
of intrinsic absorber, other magenta and blue components represent hya forest, red dots and line are observed Keck data, and black continuous line 
is the final model. 
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Figure 3: Enlarged part of Fig. 2 where the Lya region of the Al component is fitted with three different covering factors for Hi absorber. Solid- 
black line is for covering factor equal 0.36, dashed-blue - 0.61, and dotted-red - 1, i.e. complete covering. Again data from Keck are marked by red 
points and line. 
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In Fig. 2 we show that the covering factor implied for CIV component can provide good fit to the Hi component. 
Thus both covering factor and the kinematic width are the same for the two elements, and the absorption is likely 
coming from a single cloud. The fit is actually somewhat better for low covering factor although the results are not 
firm since we see well only a fraction of the line due to overlapping narrower components from Ly a forest (see Fig. 3). 

Although the errors are huge, the values of ratio are always much higher than unity, which indicates either over- 
abundance of carbon, or peculiar photoionisation model. For farther analyses we refer to the value of that ratio equal 
20, as found in Paper 1 . 



4. Photoionisation models 

In order to explain high ratio of Nav to the Nm in the absorbing medium of the HS 1603+3820 we perform 
photoionisation calculations using two codes: cloudy (Ferland et al., 1998) version 08 and titan (Dumont et al., 
2000). 

We model the broad band spectrum of HS 1603+3820 shown in Fig. 1 using two power law components, both 
with an exponential cut-off at low and high energies: 

F E = E~ a exp{-E I E max )exp{-E mi „l E), (1) 

where spectral index a = F - 1, as usual. The first power law represents the optical/UV spectrum: auv = 0.36, 
E uv = 1.36 x 10~ 4 keV, with E^ ax being first free parameter of the model, and the second power law represents the 
X-ray emission (a x =0.91, E x min = E™ x , E* ax = 100 keV). 

The important parameter is UV power-law high energy cut-off, since it regulates amount of soft photons affecting 
the absorber. For the solid line in the Fig. 1, which represents the best shape of modelled spectrum, the value E^ x = 
40 eV. But this value is not well constrained since we don't have enough far-UV observational points to determine 
overall curvature of the quasar spectrum, and since we do not consider any spectral variability. Alternatively, the 
transition between power-laws can be sharp with £^J V = 10 eV. Therefore, we adopt the value of E^ aK as a free 
parameter of the model. 

Other parameters needed for photoionisation calculations are set in a standard way. We use the ionisation param- 
eter defined as: 

where L is the quasar bolometric luminosity, no is the hydrogen number density at the cloud illuminated surface, 
and R is the distance of the intrinsic absorber from an UV/X-ray source, i.e. quasar nucleus. In all computations 
presented below, we calculate grid of models for different no and R, fixing the total quasar bolometric luminosity 
at L — 10 46 erg s _1 . For photoionisation modelling it is no important what value of L we assume, since it can be 
modified by moving cloud closer or farther from the nucleus. 



4.1. Modelling with CLOUDY 

All photoionisation models computed by cloudy assume transmission of UV/X-ray flux through a constant density 
slab, which in such case is equal to surface hydrogen density number, no- We have calculated large grid of models for 
different high energy cut-offs in UV band, E^ x , ranging from 3 eV up to 42 eV, and different ionisation parameters, 
£ from 10~ 2 up to 2.5 x 10 4 . 

The calculations were repeated for several distances, R, from the quasar nucleus, equal: log (R) = 16, 16.7, 17, 
18, and 20 [cm]. The choice of distances was made to cover all expected locations of the absorber, either close to the 
nucleus, at R — 0.00324 pc, or far away from the central black hole, at R — 32.4 pc. Since bolometric luminosity is 
the same for all models, the density of the slab is calculated from Eq. (2), and spans range from 10 2 up to 10 14 cm" 3 . 

Since column densities of ions can vary rapidly approaching the hydrogen ionisation front, located typically at 
log(A^" ) ~ 20.6 [cm" 2 ] (see Fig. 1 in Korista et al., 2008), we consider two cases of clouds. One optically thin, 
with total hydrogen column density equal to log(A^" ) = 19 [cm -2 ], and one optically thick, with column density 
log(iVjf ) = 22 [cm" 2 ]. 

In Fig. 4 we present the ratios of ionic column densities of Civ and Hi for different ionisation parameters and UV 
high energy cut-offs, for clouds with total column density log(A^" ) = 19 [crrT 2 ]. Uppermost panel shows results for 
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Figure 4: Contour plots for the ratio of Civ to Hi column densities from cloudy for clouds with total column density log( A'^" ) = 19 [cm 2 ]. 
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Figure 5: Contour plots for the ratio of Civ to Hi column densities from cloudy for clouds with total column density logt A^' ) = 22 [cm 2 ]. 
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absorber located farthest away from a nucleus at log (R) = 20 [cm]. Lower panels represent absorbers closer to the 
quasar centre, and distances are marked in their upper-left corners. The value of ionisation parameter is shown on left 
vertical axis, the E^ x on the horizontal axis, while hydrogen density number, equal «o, on right vertical axis. Colour 
contours depict the values of ratio marked in upper-right corners of each panel. 

Clouds located farther from the UV/X-ray source have low densities, and the ratio of Na v IN Hi never exceeds 1 1 . 
The highest value of Na v /Nm is achieved for absorber located at log (R) = 16.7 [cm] and reaches 18.34. 

The same grid of models but for log(A^" ) = 22 [cm -2 ] is presented in Fig. 5. Only for the dense cloud, «o = 
10 12 cm" 3 , located close to the nucleus, log (R) = 16.7 [cm], the ratio of Na v /Nm reaches 20, which agrees with our 
observational result presented in Tab. 1, first row, so we accept this model for farther considerations in Sec. 5 below. 

Finally, in Fig. 6, we present the set of clouds similar like in Fig. 4 computed with twice higher abundance of 
carbon. We see, that the distance for which the ratio of Nav IN Hi is above 20 appeared for log (R) =17 [cm], so 
the distance increased by a factor of 2. For carbon abundance five times higher than solar in each distance there is 
a possibility to achieve an absorber with Nav INni above 20, and we are not able to find one photoionisation model 
which describes observations, therefore for further comparison we use only clouds with solar metallicity. 

4.2. Modelling with TITAN 

One of the key issues in understanding of the nature of the absorbing medium is its dumpiness. A clumpy medium 
can spontaneously form out of a continuous wind, if the condition for thermal instability is satisfied. An irradiated 
medium may indeed be a subject of thermal instability, as pointed out by Krolik et al. (1981), so a cold dense cloud 
can coexist with a hot rare medium under constant pressure and the size of the cloud is roughly determined by the 
Field length (Field, 1965; Begelman and McKee, 1990; Torricelli-Ciamponi and Courvoisier, 1998; Rozahska, 1999). 
It is therefore of interest to check whether the observed systems of clouds are consistent with the requirement of the 
thermal instability. 

For that purpose, we consider the absorption by a cloud under constant pressure using the code titan, developed by 
Dumont et al. (2000). The code is designed to calculate the transmission of radiation through the cloud of considerable 
optical depth, since both the continuum and lines are calculated using the full radiative transfer instead of escape 
probability formalism. The density profile within the cloud is calculated self-consistently from the requirement of the 
total (i.e. radiation + gas) pressure being constant across the cloud (Rozahska et al., 2006). Since such numerical 
computations are time-consuming, we do not calculate an extended grid of models, but we aim to reproduce the 
observed Civ to Hi ratio of HS 1603+3820 . 

We consider a range of the total column densities from 10 21 up to 10 23 cm" 2 . Since we compute transfer through 
the cloud, we are able to make a plot of the ionisation parameter S introduced by Krolik et al. (1981): 

rn r ion * 'rati 



which shows best the thermal instability region. For each of the model computed by titan we fix cut-off energy in 
UV/EUV band on the value E^ aK = 40 eV. The free parameters of our model are the ionisation parameter and 
density number no, both at the cloud surface, and the total column density N'£'. The method is illustrated in Fig. 7 
where we show example of the temperature T vs. S relation for several models. For the purpose of this paper we have 
calculated several clouds with following sets of parameters: cloud I: no = 10 I0 cm" 3 , £, = 10 6 , = 10 23 3 cm" 2 ; 
cloud II: n = 10 10 cm" 3 , £ = 10 4 , N™ = 10 22 7 cm" 2 ; cloud III: n = 10 s cm" 3 , f = 10 6 , N™ = 10 22 cm" 2 ; cloud 
IV: n = 10 9 cm" 3 , £ = 10 4 , iVg' = 10~ 2 cm" 2 ; and cloud V: n = 10 10 cm" 3 ,f = 5x 10 3 , N'°' = 10 22 cm" 2 . 

Even for the soft spectrum, in case of HS 1603+382, the instability is still present, contrarily to the previous 
expectations (Krolik et al., 1981), but it appears at much lower optical depth since the Compton temperature of the 
incident radiation flux is considerably lower. Such a lower optical depth is in agreement with relatively low column 
densities of specific ions in comparison with those detected in a number of Seyfert 1 galaxies (Rozahska et al., 2006, 
and references therein). 

In Fig. 8, we present total column density of the cloud versus column density of Civ in the left hand panel, and the 
ratio of Nav INni in the right panel for several models. The column density of Civ is the monotonic function of total 
column density, but ratio tends to have maximum. The best representation of the observed Civ content was obtained 
for no = 10 10 , £ = 10 4 , and = 2 x 10 22 cm" 2 - marked as a red point in Fig. 8. For this model we achieved 
log(Nciv ) = 14.71 [cm" 2 ] and Na v /Nm = 19.70, which is very similar to the values observed in system Al. 
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Figure 6: Contour plots for the ratio of Civ to Hi column densities from cloudy for clouds with total column density \og(N'£' ) - 19 [cm -2 ], and 
with larger carbon abundance by a factor of 2. 
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Figure 7: The stability curve, i.e. temperature versus ionisation parameter H for constant pressure models computed by titan code. Different clouds 
are computed for the same ionisation parameter if = 10 s (see Eq. (2)) at the cloud surface, and the same total column density N'£' = 10 22 cnT 2 . 
Different hydrogen number densities on the cloud surfaces are presented by: long-dashed line - no = 10 10 , dotted line - no = 10 9 , solid line - 
no = 10 s , and dashed-dotted line - no = 10 6 cnT 3 . 




log(N CiV ) log(N clv /N HI ) 



Figure 8: Total column density versus column density of Civ - left panel, ratio of Na v INm - right panel, for different constant pressure clouds. 
Green solid line represents models for no = 10 10 cnT 3 and ^ = 10 6 , red line represents models for no = 10 10 cm -3 and £ = 10 4 , magenta line for 
no = 10 s cm~ 3 and f = 10 6 , blue line for no = 10 9 cnT 3 and if = 10 4 , and black line for no = 10 10 cm" 3 and if = 5 X 10 3 . Each computed cloud is 
represented by point on those panels. 



Note, that for constant pressure clouds, there are models presented by magenta and blue points with, «o = 10 and 
10 9 cirT 3 respectively, with ratio of Na v /Nm approaching 20, but at the same time column density of Civ is huge i.e. 
log(7Yc !V ) = 16 - 17 [cm -2 ] for those models, and exceeds the value appropriate for the Al. 

5. Distance to the Al absorber 

In the aim to estimate the distance to the Al, we have to know bolometric luminosity of the quasar. Then, the 
photoionisation modelling presented above together with the observational constrains allows us to conclude that the 
best representation of system Al is cloud with £ = 7 and «o = 10 12 from cloudy and with £ = 10 4 and «o = 10 10 
from titan. The huge difference in ionisation parameter on the cloud surface is because a cloud computed by titan 
is stratified and contains both hot and cold phases (see Fig. 7), and a cloud computed by cloudy contains only cold 
phase of the absorber. Finally, after manipulating of equation (2) we derive the distance. 

The most important is thus to estimate bolometric luminosity of the considered quasar. Misawa et al. (2005) gave 
the value of bolometric luminosity based on the prescription of Narayanan et al. (2004), where: 

Lboi ~ 4.4 AL A at X = 1450A . (4) 

For HS 1603+3820 Misawa et al. (2005) have obtained L bol = 2.5 x 10 48 ergs s" 1 . We did similar evaluation, 
besides that we have adopted the value of f° bs = 550 yuJy at A — 1450 A, as reported for this quasar by Scott et al. 
(2000) (see their Tab. 3). Assuming standard cosmological parameters: Hq = 72 kms^'Mpc -1 , Qm = 0.3, and 



11 



Table 2: Distance derivations based on cloudy and titan modelling for the Al in HS 1603+3820, for solar metallicity. Middle rows contain results 
based on bolometric luminosity calculated in this paper from Eq. (4) using flux reported by Scott et al. (2000). Bottom rows - the one calculated 
by Misawa et al. (2005). Ionisation parameter, f on the surface of the cloud drops significantly inside the cloud for the titan model, and remains 
constant for cloudy model. See text for explanation. 



Par. 


Unit 


CLOUDY 


TITAN 




[erg s~' cirT 1 ] 


7 


10 4 




[cm- 3 ] 


10 12 


10 10 


pUV 

max 


[eV] 


10 


40 


log(JVav ) 


[cm~ 2 ] 


14.97 


14.71 


Nar INm 




20.37 


19.70 


Lboi = 7.7 


xl0 4/ [ergs" 1 ] 


This paper 




log (R) 


[cm] 


17.52 


16.94 


R 


[pc] 


0.106 


0.028 


Lboi = 2.5 


xl0 4S [ergs- 1 ] 


Misawa 05 




log (R) 


[cm] 


17.78 


17.20 


R 


[pc] 


0.192 


0.051 



£2a = 0.7 for luminosity distance, and using Eq. (4), we have computed that L bo i = 7.7 x 10 erg s . This value 
is lower from that obtained by Misawa et al. (2005), by a factor of 3.25. Note, that such derivation of bolometric 
luminosity has the factor of two uncertainty. 

Taking into account both values of bolometric luminosities, we derive the distance to the Al for models computed 
by cloudy and titan. 

Regardless of the value of L ao i, in each case the Al absorber is located very close to the quasar nucleus within 
R = 0.1 - 0.2 pc. Our results differ from these obtained by Misawa et al. (2005, 2007) using variability method 
Hamann et al. (1997). Assuming that variability is caused by the change of ionisation, they found for At ~ 0.36 yr 
the location of system A within r < 6 kpc, and for 1.2 yr the location of the mini-BAL in system A within r < 8 kpc 
(Misawa et al., 2005). 

Nevertheless, our results are in agreement with dynamical model of Murray and Chiang (1995) adopted for 
HS 1603+3820 by Misawa et al. (2005) (see Sec. 6.2). Using the relation for the radius of the gas parcel relative 
to its launch radius, they have obtained constraint of r < 0.2 pc for system A. All distance derivations are given in 
Table 2. 

Our results are less restrictive if the lower value of the Civ to Hi ratio is used, or metallicity is allowed to be higher. 
The increase in the metallicity by a factor of two moves the upper limit for the distance also by a factor of 2. However, 
if we adopt lower cloud covering factor and, in consequence lower value of the Nav INm ratio, the upper limit for 
the cloud distance moves to the 30 pc. The lower limit to the distance remains practically unchanged. 

6. Connection with an accretion disk and BLR 

The distance to the Al absorber of HS 1603+3820 derived by us suggests that the absorbing material may be 
connected with the wind from an accretion disk atmosphere. Upper atmospheric layers can be quite dense up to 
riH ~ 10 14 cm 3 , depending on the distance from the central black hole as pointed by Hryniewicz (2011). Such 
density can be calculated by solving an accretion disk vertical structure parametrised by the black hole mass, its 
spin and the accretion rate (Rdzariska et al., 1999). We can properly derive density at t ~ 2/3, since we assume 
diffusion approximation for radiation when solving the disk vertical structure. This approximation is valid for the disk 
interior, but does not solve properly the structure of outer atmosphere at t <k 1. However, the density at r = 2/3 is 
representative as an initial density in the disk wind, which supplies matter to the BLR or other intrinsic absorber. 

Radial profile of the density atmosphere at the optical depth r ~ 2/3, in an accretion disk around non-rotating 
black hole in the centre of weak line quasar SDSS J0945533. 99+100950 is presented in Hryniewicz (201 1). Authors 
have shown visible bump in disk's density distribution caused by the drop of the opacity in the region where the 
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Figure 9: The value of number density on the disk surface, calculated from an accretion disk vertical structure (see Sec. 5 for explanation), versus 
the distance from the quasar centre. Two radial profiles (blue and black) are calculated for two viscous heating prescription marked in the left upper 
corner. Points represent values achieved for both parameters from photoionisation modelling. Clouds computed with cloudy (c), are marked by 
circles, and those computed with titan (t) - by triangles. Red points represent clouds with bolometric luminosity calculated by us, while green 
points - by Misawa et al. (2005). 

temperature is too low for significant contribution of electron scattering and free-free processes (Hure et al., 1994) 
but still too high for opacity due to the dust, so only molecules are the source of opacity. Sudden drop of the volume 
hydrogen density from 10 14 to 10 9 cirT 3 occurs at log (R) ~ 18 [cm] (see Fig. 5 in Hryniewicz (201 1)), and coincides 
perfectly with the distance to BLR for this object. 

We have made the same calculations for the mass of the black hole appropriate for HS 1603+3820 , derived from 
bolometric luminosity L/, / = 7.7 x 10 47 erg s _1 , assuming mass accretion rate equal Eddington, black hole mass is 
Mbh - 5.26 x 1O 9 M . Furthermore, we have estimated black hole mass in this quasar using Civ line (Vestergaard and 
Peterson, 2006), which gives similar value Mbh — 5.4 x 1O 9 M , with uncertainty of one order of magnitude. 

Taking into account bolometric luminosity derived by us, we consider an accretion disk around black hole of the 
mass 5.26 x 10 9 M o , and we assume Eddington accretion rate for the gas. For such disk we have computed radial 
profile of the density at r = 2/3, assuming two different viscous heating prescriptions. First, we assume that 
viscous heating is proportional only to the gas pressure, so we do not expect any radiation pressure instabilities. 
Second, we assume that viscous stress is proportional to the total pressure. 

Fig. 9 presents radial profiles of for two different viscosity prescriptions. The drop in the density of the disk 
atmosphere at the radius log(^) ~ 17 [cm] is related to the partial ionisation of the hydrogen in the disk atmosphere 
and a very strong inversion in the vertical density profile due to radiation pressure. This partial ionisation effect is not 
related to the overall disk instability since the ionisation instability operates when there is partial ionisation of the disk 
interior (i.e. further out). What is more, if we change the viscosity prescription from aP rot to aP gas , the plot practically 
does not change, although the integrated disk surface density, 2, changes more than by an order of magnitude. Such a 
strong density inversion in the disk atmosphere can lead to local instabilities and an outflow. Additionally, all points 
from our modelling are marked on this diagram. They are in good agreement with values predicted for outer layers of 
an accretion disk atmosphere around black hole with mass estimated by us from observational constrains. 

Furthermore, following equation (3) in Pian et al. (2005), we can calculate the distance to the BLR in HS 1603+3820, 
from the continuum luminosity at 1350 A. Estimation is rather crude, since absorption can be present at this wave- 
length, but we have got the value of log(RBut) = 18.75 [cm]. Together with our distance derivation to the Al 
absorber, we conclude that intrinsic absorption is located in radial direction closer to the nucleus than a possible BLR 
in HS 1603+3820. 



13 



7. Conclusions 

We have presented photoionisation modelling of intrinsic absorber in quasar HS 1603+3820. Using observations 
from optical/UV/X-ray energy bands, we have constructed the shape of continuum, which interacts with absorbers on 
the way from the nucleus. Double power-law with exponential cut-offs and related observed optical to X-ray slope 
a M = 1 .49, was adopted as an input to photoionisation modelling of the first cloud in system A. 

For several lines of system Al, reported in Paper I, we have found very high ratio of column densities of Civ to Hi 
. The value of this ratio in some cases exceeds 20. Our photoionisation computations done by cloudy and titan show 
that this value, together with intrinsic SED of HS 1603+3820, allowed us to choose a single solution which explains 
all observational constrains. 

From modelling done by titan, we found that total column density of absorbing system is 2 x 10 22 cm -2 . This 
value is quite close to the maximum column density in a constant pressure medium due to the thermal instability. It is 
one order of magnitude lower that the values derived for Seyfert galaxies (Rozariska et al., 2008), since the intrinsic 
spectrum of HS 1603+3820 is much softer, with a very weak X-ray tail. This coincidence supports the view that 
intrinsic absorbers do form in the wind as a result of thermal instability and the constant pressure model is better 
suited to interpret the data than the constant density one. 

The hydrogen number density derived from modelling is rather high, of the order of n = 10 10 , 10 12 cm -3 , which is 
consistent with previous density estimations of UV intrinsic absorbers in NLS1 galaxies (see Leighly (2004)). On the 
other hand, our result is inconsistent with number density found by Misawa et al. (2005) using variability method. 

Such high value of number density gives the location of absorbing gas quite close to the nucleus within R = 
0.1 - 0.2 pc. This location supports scenario that UV/X-ray outflows are in the form of winds from accretion disk 
atmospheres with similar number densities. Furthermore, we have shown that the Al absorber is located closer to the 
nucleus than the BLR for a given bolometric luminosity of HS 1603+3820. 

The results, however, are very sensitive to the measured Civ to Hi ratio, which is sensitive to the determination of 
the covering factor, very difficult for hydrogen line. For the same covering factor for Hi and Civ of 0.36 in the Keck 
2001 data, the ratio is much lower than measured in Paper I and the upper limit for the distance is moves to 30 pc. 
Also assumption of higher metallicity increases the upper limit for the absorber distance. 

Close location of the absorber system is in contradiction with Misawa et al. (2005, 2007) who found that sys- 
tem A is located within r — 6 kpc, but it agrees with dynamical model of (Murray and Chiang, 1995) adopted for 
HS 1603+3820 by (Misawa et al., 2005) (see Sec. 6.2). Using the relation for the radius of the gas parcel relative to 
its launch radius, they have obtained constraint of r < 0.2 pc for system A. 

The outflow velocities of system A are rather high, up to few thousands km s , which is consistent with the winds 
detected in X-rays in several NLS1. Therefore, HS 1603 +3820 can be considered as a high redshift analog of NLS1 
galaxies. 
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